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We studied solar-like oscillations in 115 red giants in the three open clusters 
NGC 6791, NGC 6811, and NGC 6819, based on photometric data covering 
more than 19 months with NASA's Kepler space telescope. We present the 
asteroseismic diagrams of the asymptotic parameters Si>02, ^oi and e, which 
show clear correlation with fundamental stellar parameters such as mass and 
radius. When the stellar populations from the clusters are compared, we see 
evidence for a difference in mass of the red giant branch stars, and possibly a 
difference in structure of the red clump stars, from our measurements of the 
small separations <5z/o2 and Svqi- Ensemble echelle diagrams and upper limits to 
the linewidths of I = modes as a function of Av of the clusters NGC 6791 
and NGC 6819 are also shown, together with the correlation between the i = 
ridge width and the T eS of the stars. Lastly, we distinguish between red giant 
branch and red clump stars through the measurement of the period spacing of 
mixed dipole modes in 53 stars among all the three clusters to verify the stellar 
classification from the color-magnitude diagram. These seismic results also allow 
us to identify a number of special cases, including evolved blue stragglers and 
binaries, as well as stars in late He-core burning phases, which can be potentially 
interesting targets for detailed theoretical modeling. 

Subject headings: open clusters and associations: individual (NGC 6791, NGC 6811, 
NGC 6819) - stars: oscillations - stars: evolution - techniques: photometric 



Introduction 



Many recent studies of solar-like oscillations of red giant stars have been focused on large 
ensembles of stars, made possible by the flood of high quality photometric data provided 



(e.g. 


De Ridder et al. 


2009 


Kallinger et al. 


2010b 




Vlosser 


Borucki et al. 


2010 


Koch et al. 


2010 


Bedding et al. 


2010 


)• 



Particular attention has been given to the three open clusters NGC 6791, NGC 6811, and 



NGC 6819 in the Kepler field (Stello et al. 2010 Basu et al. 2011 Hekker et al. 2011 Stello 



et al. 2011a[b ; Miglio et al. 2012), due to the well-known advantage of cluster stars sharing 
common properties, which allows for more stringent investigations of stellar evolution theory. 

Among the highlights in recent results relevant for our study are the measurements of 
the small frequency separations 5^02, Suqi and of the dimensionless term e, their correlation 



Bedding et al. 2010 


Huber et al. 


2010 


Kallinger 



et al. 2012, for previous results on field red giant stars) and their dependence on stellar 



mass ( Montalban et al7||2010 Kallinger et al.||2012 ). Also, the results on the ensemble echelle 
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diagrams have allowed for the investigation of ensemble properties of the modes, including the 
measurement of the mean small spacing 5i/ 03 dBedding et al.|2010[ |Huber et al.|[20lo| and the 
linewidths of the dipole modes and their correlation to fundamental stellar properties (e.g. 



sec 



Chaplin et al. 2009 Baudin et al. 2011 Appourchaux et al. 2012 Belkacem et al. 2012), 



which are important for the comprehension of the physics responsible for the excitation and 
damping of solar-like oscillations. Finally, the period spacing analysis for the investigation 



of the evolutionary stage of red giants (Bedding et al. 2011 Mosser et al. 2011a) now allows 



us to distinguish between He-core burning red giants and those only burning hydrogen in a 
shell. 

Here, we study 115 red giants belonging to the above mentioned clusters, continuously 



observed for 19 months by the NASA Kepler Mission (see Jenkins et al. 2010b Gilliland et 



al. 2010 Brown et al. |20lT , for details on the data pipeline and acquisition, for a general 
introduction to the asteroseismic program, and for a description of the Kepler Input Catalog, 



respectively). Our study is made along the same lines as those described by Huber et al. 
(2010), who analyzed the first 4.5 months of Kepler observations of field red giants. 



In 



contrast to 



Huber et al. (2010), our cluster red giants have the great advantage of providing 



more homogeneous samples because age, metallicity and mass (see Basu et al. 2011 Miglio 



et al.|2012[ ), are about the same. In particular, NGC 6791 is a very old, (~ 8.3 Gyr, |Brogaard 



et al. 2012), and metal rich ([Fe/H] = 0.29 ± 0.03 (random) ±0.07 (systematic), Brogaard 



open cluster, with average masses M RGB = 1.2O±O.O1M ( |Basu et al.pOlT] ) and 
M RC = 1.15 ± 0.03 M Q (|Miglio et gpOjjb for red giant branch (RGB) and red clump (RC) 



stars respectively (see also |Brogaard et aL 2012, for recent results from eclipsing binaries). 
NGC 6819 is a middle aged (2-2.4 Gyr, |Basu et~aL"]2011fr open cluster, with solar metallicity 
([Fe/H] = 0.09 ± 0.03, |Bragaglia et~aT||200l| ), and average masses M RGB = 1.68 ± O.O3M 
and M RC = 1.65 ± 0.04 for RGB and RC stars respectively. The third open cluster, NGC 



6811, is characterized by a young (0.7 ± 0.1 Gyr, Glushkova et al. 1999) and possibly solar 



metallicity star population (suggested by two independent spectroscopic investigations by 
Bruntt et al., in prep, and Molenda-Zakowicz et al., in prep.), where a small number of RC 
stars has been observed, showing an average mass M RGB = 2.35 ± 0.04 M (Stello et al. 
2011a[b| . The temperature estimates for both NGC 6791 and NGC 6819 were derived by 



Hekker et al. (2011). In particular, they used color-temperature calibrations by Ramirez & 



Melendez (2005) and JHK photometry from the 2 MASS catalog (Skrutskie et al. 2006), 



which is available for all the stars of the sample. V magnitudes are taken from Stetson 



et al. (2003) for NGC 6791 and from Hole et al. (2009), in order to derive temperatures 
based on the (V — K) color. The adopted reddenings are E(B — V) = 0.16 ± 0.02 for 



NGC 6791 (Brogaard et al. 2011) and E(B - V) = 0.15 for NGC 6819 (Bragaglia et al. 



2001). Lastly, Basu et al. (2011) estimated DMs for both NGC 6791 and NGC 6819 by 
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adopting the extinction Ay = 3.1E(B - 
(m — M) = 11.85 ± 0.05, respectively. 



V), which yielded to (m - M) = 13.11 ± 0.06 and 



After briefly introducing the parameters involved in our study in Section |2j we describe 
the code developed for this work in Section [3j which concerns the analysis of the average 
p-mode structure in the power spectrum. In Section [4] we show the resulting asteroseismic 
ensemble diagrams and the linewidths of radial modes as a function of fundamental stellar 
parameters, while Section [5] presents the analysis of the period spacing of mixed dipole modes 
following the approach of Bedding et al. (2011). Finally, we conclude in Section [6] 



2. Asymptotic parameters 



Before proceeding with the description of the data analysis, it is helpful to introduce 
the physical quantities that we study in this work. As known from the asymptotic theory 
of solar-like oscillations, acoustic standing waves (also known as p modes) with low angular 
degrees, £, and high radial orders, n, show regular frequency spacings according to the 



^Vandakurov 


1968 


Tassoul 


1980 


Gough 


1986) 



Here 



v nl ~ Av 



Au 



n + - + e ) - 5v e . 



dr 



(2) 



'o c(r) 

is known as the large frequency separation, which scales roughly as the square root of the 
mean stellar density (Ulrich 1986), while e is a phase shift sensitive to the properties of the 



near surface layers of the star ( Christensen-Dalsgaard & Perez Hernandez 1992). The term 
5is £ is known as the small frequency separation and for £ = 1, 2 and 3 it is defined as 



8u 01 



5^02 
1 

~~ 2 
1 

" 2 



(^n,0 + ^n+1,0 
(^n,0 + ^n+1,0 



) - Vn,l, 



*V3- 



(3) 
(4) 
(5) 



The small frequency separations are related to the sound speed gradient in the stellar core, 
hence to the mean molecular weight, which increases as the star evolves. Mixed modes occur 
as the frequencies of the g modes in the core and the p modes in the envelope become similar 
during the subgiant and red giant phase. As the star evolves, its mixed modes will undergo 
avoided crossings causing so called mode bumping, which broadens the ridges in the echelle 
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diagram (e.g. see Stello 2011). As argued by Bedding et al. (2010) and Huber et al. (2010) 



the small spacing 5u 03 is therefore preferred in red giants over the more conventional 81/13. 
We note that mode bumping mostly affects the dipole modes as they penetrate deeper into 



the star, and hence couple more strongly to the g modes in the core (Dupret et al. 2009). 
The dipole modes are therefore sensitive to the core properties of the star, which allows us 
to determine which red giants burn helium or not (see Section [5]). 

Figure [T] shows a typical power spectrum of a low luminosity RGB star, KIC 2436593, 
observed in NGC 6791. The mode identification for some of the modes is shown, together 
with the indication of regions containing mixed i = 1 modes, represented by the gray-shaded 
strips. The large separation Av and the small spacing 5i/ 2 are indicated as well. The inset 
shows a zoom-in of one of the gray-shaded strips, where the observed period spacing of the 
dipole modes, AP Q b s is marked. 



3. Observations and data analysis 



The photometric time series of the 115 red giants used in this work were obtained in 
Kepler 's long cadence mode (At ~ 30 min, Jenkins et al. 2010a) between 2009 May 13 



and 2010 December 22. This corresponds to the observing quarters 1—7, providing a total of 



almost 18,000 data points per star (see Garcia et al.||2011 for details on the detrending of the 



data). We followed the approach described by Stello et al. (2011b) for merging the quarters, 
and we discarded the stars that they classified as seismic non-members in their study. We also 
note that, according to their classification based on the color- magnitude diagram (CMD), 
the cluster stars in NGC 6811 are all He-burning stars, with one star appearing to be in a 
late He-core burning phase towards the asymptotic giant branch (AGB). 

For the present study we developed the Asymptotic Analysis of Red Giants (AARG) 
code, with the purpose of deriving asymptotic parameters for p modes and observed period 
spacings for mixed modes in red giant stars. AARG performs a multi-step analysis in a semi- 
interactive way, allowing the user to follow the results at each step and make any necessary 
corrections. We calculated background-corrected power spectra and measured Av using the 



SYD pipeline (Huber et al. 2009). As a check we compared Av values with those derived 



using other methods (Mosser & Appourchaux 2009 Mathur et al. 2010 Hekker et al. 2010 



Kallinger et al.||2010a ) and found good agreement. We focus first on the analysis of p modes, 
which represents the main part of the work, leaving the discussion of period spacings to 
Section HI 



The analysis of p modes, performed for each star, was done in three steps: (i) collapse the 
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echelle diagram using the measured Au; (ii) identify the centroids z/ , u 1 , v 2 of the £ = 0, 1, 2 
ridges by fitting three Lorentzian profiles to the collapsed echelle diagram, which gives the 



small spacings 5v 2 and 8voi, and e (see the next paragraph and Huber et al. 2010); and 
(iii) simulate 500 power spectra by perturbing the observed power spectrum of the star 
according to a x 2 statistics with 2-degrees of freedom, perform the first two steps of the 
analysis for each simulation in order to derive a new set of asymptotic parameters, and 
evaluate their uncertainties by computing a robust rms of the results. Figure [2] shows an 
example of a collapsed echelle diagram obtained with the AARG code. The centroids of the 
ridges £ = 0,1, and 2 are marked by dotted lines, while the Lorentzian profiles used to fit 
the different ridges are shown with thick solid lines (red, blue, and green, respectively). 

We now describe step (ii) in slightly more detail. AARG requires an initial guess of e 
for each star, given by manually marking the position of the £ = peak in the collapsed 
echelle diagram. This is followed by a Lorentzian fit to obtain the centroid, vq. According 



to the values shown by Huber et al. (2010) and White et al. (2011b), we expected to have 



0.5 < e < 1.5 for red giants with Au < 15 /iHz. We therefore either added or subtracted 1 to 
the measured e to ensure it would fall within this range (see also |Mosser et af|2011b| jWhite 



et al. 2011a, for further discussion of the position of the centroids and the e diagrams). 



Next, AARG makes a first guess for the centroid v 2 of the £ = 2 peak by adopting an 
empirical relation u — v% = 5u 02 = cAu, where c is small. Although the relation has a slight 
mass dependence (see Montalban et al. 2010), using a fixed value of c = 0.123 offered a 



reliable first guess for every star in our sample (note that our value is very close to the one 



measured by Bedding et al. 2010). As for u , a Lorentzian fit centered on the first guess for 
the £ = 2 ridge position provides the final value of the centroid u 2 , and hence also Su 02 . The 
search for the £ = 1 peak is performed automatically by finding the maximum in the regions 
of the collapsed echelle diagram laying outside the £ = 0,2 peaks. A third Lorentzian fit is 
then performed, providing the centroid v\, which gives Svqi = Vq + Au/2 — v\ according to 
the convention by Bedding et al. (2010). For a few stars (~ 10) our method did not perform 
well. This was mainly caused by partly overlapping £ = 0, 2 peaks (especially in NGC 6811 
because of the higher mass of its stars) and strongly affected £ = 1 peaks due to mixed 
modes. These stars were manually analyzed afterwards. We could successfully derive the 
asymptotic parameters of p modes for a total of 115 stars: 60 for NGC 6791, 5 for NGC 6811, 
and 50 for NGC 6819. 
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350 pi 




Frequency (//.Hz) 



Fig. 1. — Power spectrum of KIC 2436593, a typical low luminosity RGB star belonging 
to NGC 6791. Mode identification for some of the peaks is shown. Shaded regions in gray 
indicate mixed i = 1 modes. Av and 5^02 are also marked. The inset shows the detail of 
one of the i = 1 shaded regions, where an indication of the observed period spacing, AP Q b s , 
is shown. 
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Fig. 2. — Collapsed echelle diagram of KIC 2436593. The identification of the ridges £ = 
0, 1, and 2 is shown, together with their Lorentzian fits (red, blue, and green solid lines, 
respectively). 




-8- 



Results 



4.1. e diagram 

The e term of Equation was shown to be highly correlated with Av for red giant 



stars (Huber et al. 2010 Mosser et al. 2011b). The e diagram is shown in Figure 3 for the 



clusters NGC 6791, NGC 6811, and NGC 6819, where 1-a error bars were derived by means 
of Equation Q. We note that the RC stars (identified from the CMDs by Stello et al. 



(2011b) but adjusted for a few stars based on our analysis of the period spacing, presented 
in Section [5]) form distinct groups with slightly lower e than the RGB stars at Av ~ 3.7 //Hz 
for NGC 6791, Av ~ 8 //Hz for NGC 6811, and Av ~ 4.8 //Hz for NGC 6819. In particular, 
we measured a weighted average of e for clump stars and RGB stars in the same Av range 
of RC stars, and found them to be (e RC ) = 0.829 ± 0.031 and (e RGB ) = 0.915 ± 0.039 for 
NGC 6791, and (e RC ) = 0.970 ± 0.018 and (e RGB ) = 1.015 ± 0.017 for NGC 6819. In both 
cases, appears to be significantly different from (e RGB ). This is in good agreement 



with Bedding et al. (2011) and Kallinger et al. (2012). Although a lower mass of RC stars 



can result in lower Av, and hence lower e, one should note that this effect alone cannot 



explain the observed difference in e between RC and RGB stars of similar Av (see Miglio 



et al. 2012, for a detailed study about the mass difference between the RGB and RC stars). 
The difference in evolutionary state also needs to be taken into account to fully explain the 



observed difference in phase shift (Kallinger et al. 2012). 



A least-squares fit to the RGB stars of the clusters was computed, using the log-relation 



e = A + 5 log Av. 



(6) 



adopted by Mosser et al. (2011b). Since the fits computed to the RGB stars of NGC 6791 



and NGC 6819 are not significantly different, we give the result for all the RGB stars in our 
sample, providing single values for the coefficients A and B. The result is shown as a solid 
black line in Figure |3j where A = 0.601 ± 0.025 and B = 0.632 ± 0.032. The fit from |Mossef 



et al. (2011b), who used a five-month data set, is added for comparison and plotted in the 



Av range [0.6, 10] /iHz covered in that study (dashed purple line). The fit by Kallinger et al. 



(2012), based on more than 900 field red giants observed by Kepler for about 600 days, is 



almost indistinguishable from ours (dot-dashed cyan line). We also tested a power-law form 
of the e-Av relation and found the x 2 to be very similar to that derived from Equation (J6|. 



The log-relation was finally chosen to allow a direct comparison with the results by Mosser 



et al. (2011b) and Kallinger et al. (2012) 



Lastly, we note that the uncertainties on Av, and hence on e, become quite large for 
values of Av below < 2 /iHz due to the limited frequency resolution and small number of 
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orders observed. For the star with the highest Av, the large uncertainty is caused by the 
low S/N level, due to its low oscillation amplitude. 



4.2. C-D diagrams 



In the C-D diagram one plots the small spacing <5z/q2 versus the large spacing Av 



(Christensen-Dalsgaard 1984), which for MS stars enables one to discriminate stars of dif- 



ferent age and mass. A new version of the C-D diagram proposed by Mazumdar (2005) and 



by Montalban et al. (2010) for MS and RGB stars respectively, is constructed by considering 



5vqi instead of 5vq2 (see also White et al. 2011a). It has been shown that for red giants the 



C-D diagrams can not be used to investigate age (White et al. 2011b) but that it is still 



useful to discriminate mass (e.g. see Bedding et al.||2010 Huber et al.||2010 Montalban et al. 
2010 Kallinger et al.|2012 ). The results for both £z/ 2 and 5vqi are shown in Figures |4]^a) and 
(c), for the three clusters. As before, open symbols represent RC stars while filled symbols 
are RGB stars. As a first approximation, we represent the relation between 5vq2 and Av by 
the linear relation 5v$2 = ao2 + &02 Az/, which we fitted with a standard least-squares method 
to the RGB stars. The results are shown in Fig. [4] with a dashed red line for NGC 6791 and 
a dot-dashed green line for NGC 6819. Their equations are given by 



5v, 



(6791) 
02 



(0.121 ± 0.003) Av + (0.035 ± 0.012) /iHz . 



(7) 



and 



5v, 



(6819) 



> 02 - (0.114 ±0.003) Av + (0.019 ±0.012) ^Hz . (8) 

Only error bars on 5v 02 were considered for the fits, but the results obtained by including 
uncertainties on both quantities were indistinguishable from the ones presented here. The 
coefficients of the 8vq2-Av relation estimated from our fit agree within a few percent with 



those derived by Kallinger et al. (2012) for field stars. 



The typical mass for an RGB star, M RGB , is expected to be different for each cluster 
but about the same within a given cluster. For cluster RGB stars we therefore have a 



much tighter constraint on the stellar mass than for field stars. Basu et al. (2011) found 
M 679 i = 1.20 ± 0.01 M and M 68 i 9 = 1.68 ± 0.03 M as the averages for the RGB stars, 
which were derived using grids of stellar models that incorporated scaling relations for v max 



and Av. We refer to Miglio et al. (2012) for further discussion about the mass estimates 



for these stars. Following the theoretical work by Montalban et al. (2010), who showed that 



Su 2 depends on mass, we relate the difference in the slopes, 602, in Equations ^ and ^ to 
the difference in Mrgb- Assuming the linear relation 

'M RGB 



J()2 



«02 + A 



02 



Mrr 



(9) 
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Fig. 3. — e diagram for the clusters NGC 6791 (red circles), NGC 6811 (blue diamonds), 
and NGC 6819 (green squares). Open symbols represent RC stars while filled symbols are 
RGB stars. 1-a uncertainties are displayed for both quantities. A fit to the RGB stars of all 
the clusters using Equation ^ is added (solid black line), as well as the one from Mosser 



et al. (2011b) (dashed purple line) and Kallinger et al. ( |2012 ) (dot-dashed cyan line). Stars 



marked with labels and arrows represent special cases that are discussed in Section 5.2 
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Fig. 4. — (a), (c): C-D diagrams of the small spacings 5z/ 2 and 5v m for the clusters 
NGC 6791 (red circles), NGC 6811 (blue diamonds), and NGC 6819 (green squares). Open 
symbols represent RC stars while filled symbols are RGB stars. Error bars show l-a uncer- 
tainties. The linear fits to the RGB stars are shown for both NGC 6791 (dashed red line) 
and NGC 6819 (dot-dashed green line). Stars marked with labels represent special cases 
that are discussed in Section 5.2 (b), (d): modified C-D diagrams of the ratios 81^02/ Av 
and Svoi/Av with the same notation adopted for the upper panels. 
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we obtain a 2 = 0.138 ± 0.012 and f3 02 
equations (one for each cluster). 



-0.014 ± 0.008 by solving the system of two 



As done for 5u 2, we fitted a linear relation 5u m = aoi + b^iAu to the RGB stars 
of NGC 6791 and NGC 6819. The linear trends are shown in Figure g}c) with the same 
notation as Figure [5](a), and the results are 



5u, 



(6791) 



01 



and 



5u, 



(6819) 
01 



(0.021 ± 0.003) Av - (0.063 ± 0.011) /iHz 



(0.000 ± 0.003) Ais - (0.109 ± 0.012) /iHz . 



(10) 



Once again, the uncertainties are quite similar for the two clusters. Like 5u Q 2, we also see 
a mass dependence on <5i/ i for the RGB stars. But unlike <5fo2, the trend appears to go in 
the opposite direction, with higher 5i/qi for higher masses. This is in qualitative agreement 



with the theoretical results by Montalban et al. (2010), whose Figure 5(b) shows a slight 
increase in <5z/ i for increasing mass along the RGB. As for 5^02, we relate the slopes 6qi m 



Equations ( JlO[ ) and (JTTj) to Mrgb, assuming the linear relation 

'M RGB ' 



b i = a 0i + ft 



01 



(12) 



and find a 01 = -0.073 ±0.012 and /3 i = 0.044 ±0.008. We find that \/3 i\ ~ 3|/3 02 |, hence it 
appears that 5u 01 is more sensitive to mass than 8u 02 by about a factor of three. But at this 
stage we would caution overinterpretation of this result as further theoretical investigations 
are required to fully understand how 5uqi depends on the fundamental parameters and 



internal structure of red giants (Section 4.3). 



Figures |4^b) and (d) show the so-called modified C-D diagrams which plot the relative 
ratios S^/Au and Suoi/Au. The reason for considering the ratio Suqi/Au is that models 



show it to be less sensitive to surface layer effects (e.g. see White et al. 2011a) and that the 
small spacings Suqi approximatively scale with Av. Our results appear to be in agreement 



with previous results on red giants ( 


Bedding et al. 


2010 


Huber et al. 


2010 


Mosser et al. 


2011b 


Kallinger et al. 


2012 


) and with the theoretical studies by 


Montalban et al. 


(2010 


)• 



4.3. The small spacings of red clump stars 

It is interesting to compare the average small spacings for the RC stars relative to the 
RGB stars in each cluster. In the following we denote this quantity A{5u^) = {5u^p) — 
(<5^ GB ). It is evident in all four panels of Figure El but slightly more so in Figures E[b) and 
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(d), that the RC stars on average show different small spacings than RGB stars of similar 
Av. We will first discuss A(5i , 02 ). 

For NGC 6819, A(6u 02 ) = 0.112 ± 0.016 //Hz, while for NGC 6791 we have A(5u 02 ) = 



0.012 ± 0.021 /tHz. Given the relation between Suq 2 and mass for RGB stars (Section 4.2), 
one might speculate that a similar relation would exist for RC stars. However, we note that 
A(8u 02 ) for NGC 6819 is about twice as large as the difference in Su 02 between the two RGB 
populations in NGC 6791 and NGC 6819 at a similar Av. Hence, if we applied the relation 
in Equation ^ to the RC stars, we would find that the RC stars in NGC 6819 have a mass of 
about O.7M (corresponding to a mass loss of about 1 M ), in stark disagreement with the 



results by Miglio et al.| (2012), who found A(M) = —0.03 ± 0.04. Hence, there is certainly 



something else dominating the different values of A(5u 02 ) we see for the two clusters. 

Turning our attention to the other small spacing, we have A(5voi) — 0.004 ± 0.025 /iHz 
for NGC 6819 and A(6u 01 ) = -0.113 ± 0.020 /iHz for NGC 6791. Hence we see that 
A(<5z/ i)679i < A(<5z/ i)68i9; which was also the case for A(<5z/ 02 )- We recall that the mass 



dependencies of Si/ 02 and 5u m were opposite for the RGB, both in observations (Section 4.2) 



and models (Montalban et al. 2010). Hence, we would also expect an opposite trend for 
A{5vqi) (A(5z/oi)679i > A (5z/oi)68i9) • The fact that we do not observe this is further evi- 
dence that a simple relation with mass alone cannot explain the observed differences in small 
spacings between RC and RGB stars. 

A possible explanation is the internal structural changes of the stars that occur during 



the He-flash phase (Bildsten et al. 2012) between the tip of the RGB and the RC. These 
changes could be significantly different for stars of different masses (M 679 i = 1.20 ± 0.01 M 
and M 6819 = 1.68 ± 0.03 M Q , |Basu et al.| (|20TT|)), composition ([Fe/H] 6791 = 0.29 ± 0.03 



(random) ±0.07 (systematic), Brogaardet al. (2011); [Fe/H] 6819 = 0.09 ±0.03, Bragagliaet al. 



(2001)), and rotation rates (Meibom et al. 2011). Further investigation requires modeling of 



both the RGB and RC stars in these clusters. 



The dependence of 5u 01 on stellar properties was investigated by Montalban et al. (2010) 



using stellar models covering 0.7—2.3 M on the RGB and 2.5—5.0 M G in the He-core burning 
phase. They found that small values of Svqi were predominantly seen among RGB models, 
and we would therefore expect the RC stars to show larger Suqi on average, which is in 
contrary to what we observe for NGC 6791. However, we note that all the He-core burning 



models in the Montalban et al. (2010) sample were more massive than the stars in the two 
clusters considered here, and a direct comparison is therefore not possible. The physical 



cause of a lower value of 5pq\ was not firmly established by Montalban et al. (2010), but 



they argued that there was a tendency for low 8vq\ values in models where the inner turning 
point of the I = 1 modes was well inside the convective envelope, corresponding to stars 



-14- 



with deep convective envelopes. Clearly, these issues deserve further study. 



4.4. Ensemble echelle diagrams 



Following Huber et al. (2010), we computed the so-called ensemble echelle diagrams 
for both NGC 6791 and NGC 6819. When dealing with a large number of stars, ensemble 
echelle diagrams are very helpful for studying the evolution of features such as ridge width 
and position. In particular, the measurement of the average position of the £ = 3 ridge 
becomes possible also when low S/N in the power spectra does not allow one to make a clear 
detection of the corresponding peak in a single star. The results are shown in Figures [5|a) 
and (c), where the stars are numbered by increasing Av. Each row in the plot represents the 
collapsed echelle diagram of a single star using the large separation adopted in the analysis, 
and shifted in order to have the £ = ridge fall on [yjAv mod 1) = 0.3 (see also the 



discussion by Stello 2011, Sec. 2). The RC stars (red star symbol) clearly show strong 
broadening of the £ = 1 and 2 ridges. We see that even the £ = ridge appears broader for 
RC stars in both clusters. 

The result of collapsing the ensemble echelle over the entire sample of stars is shown 
in Figures [5]^b) and (d) (thick black line). Results for RC stars (red line) and RGB stars 
(blue line) are also plotted for both clusters. The presence of an £ = 3 peak becomes 
evident for NGC 6791, while for NGC 6819 a hint of £ = 3 is visible only for the RGB stars. 
For NGC 6791, the £ = 3 hump seems to arise from several stars, particularly those with 
Av < 7 /iHz, as visible from Figure |5^a). The position of the marked £ = 3 peaks of the two 
clusters, and hence their average small spacings Sv 03 , are in agreement with the results of 



Bedding et al. 


(2010 


); 


Huber et al. 


(2010 


); 


Mosser et al. 


(2011b 


) and 


Kallinger et al. 


(2012 



It is noticeable that the £ = 1,2 ridges move away from the £ = ridge as the stars 
evolve from H-shell to He-core burning red giants, a result that was already discussed by 

on the left slope 



Huber et al. (2010). We also note that the hump visible in Figure [5| 
of the £ = 1 peak (red line), is caused by only two stars having strong peaks that occur at 
[y I Av mod 1) ~ 0.7 and this is therefore not an indication of a general feature. Referring 
to the effect on 5v i discussed in Section |4.3[ we notice that the £ = 1 ridge of the RC stars 



of NGC 6791 (Figure |na)) is shifted towards the right-hand side of the diagram, i.e towards 
lower values of the small spacing, while this shift is not apparent in NGC 6819 (Figure p[c)). 
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Fig. 5. — (a), (c): ensemble echelle diagrams of the clusters NGC 6791 and NGC 6819, 
respectively, where £ = centroids were aligned by shifting the £ = ridge of each star 
to align with {yjAv mod 1) = 0.3. The number of the stars, ordered by increasing Az/, is 
shown on the left coordinate, and the corresponding Az/ is shown on the right axis. Red star 
symbols mark the clump stars identified in the clusters. Note that each row corresponds to 
the collapsed echelle of one star, normalized to unity, (b), (d): diagrams showing panels (a), 
(c) collapsed over the entire range of Az/ (thick black line) normalized to unity. Results for 
RC stars in red and RGB stars in blue are also shown. Ridge identifications and definitions 
of small separations are indicated. In both panels, the dotted lines represent the centroids 
of the £ = 0, 1, 2 and 3 ridges, while the dashed line is the position of the midpoint of two 
adjacent £ = modes. 
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Fig. 6.— FWHM of the t = ridge as a function of Av for RGB stars in NGC 6791 (filled 
red circles) and NGC 6819 (filled green squares). Each point represents the average of values 
within a subset of stars with similar Aza The error bars are the l-a uncertainties on the 
mean for each subset. Open symbols at Av ~ 3.7/iHz and Av ~ 4.8 //Hz represent the 
measurements for the subsets of RC stars. 
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Fig. 7.— FWHM of the £ = ridge plotted against T eS for the stars of NGC 6791 (red circles) 
and NGC 6819 (green squares). Also shown are measured linewidths for MS and subgiant 
field stars (blue diamonds) from Appourchaux et al. (2012). Each cluster point represents 
the same subset of stars plotted in Figure [6] The error bars are the l-a uncertainties on the 
mean for each subset. The fit to the MS and subgiant stars taken from |Appourchaux et al. 



(2012) is also shown (dot-dashed blue line). The dashed black line shows an exponential fit 
(Equation (13)) to all stars. 
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4.5. Mode linewidths 



Measuring the linewidths of p modes and studying how they correlate to the fundamental 
stellar properties has important consequences for the understanding of the damped nature of 
solar-like oscillations. In fact, the physics responsible for the damping mechanism that acts 



in the convective envelope of low mass stars is not yet fully understood (e.g. see Houdek et 



al. 


1999 


Dupret et al. 


2009 



In the present work we provide estimates of the linewidths of radial modes derived 
through the AARG code. In particular, the widths of the ridges in the collapsed echelle 
diagrams (Figures |5^b) and (d)) give a rough estimate of the mode linewidths. Figure [6] 
shows the FWHM for the I = ridge from the Lorentzian fit to the corresponding peak in 
the collapsed echelle diagram, for the RGB stars of both NGC 6791 and NGC 6819. Each 
point is the average from a subset of stars sorted in bins of Av, while the overlaid error 
bars are 1-a uncertainties on the mean for each bin. The open symbols at Av ~ 3.7/iHz 
and Av ~ 4.8 //Hz show our measurements for the RC stars. We see a clear increasing 
trend when moving to higher Av, a result that was already apparent from the analysis by 



Huber et al. (2010) of field red giants, despite of the shorter data set that was available (see 



also Kallinger et al. 2012, who obtained a similar result by using a different method). This 
increasing trend is also visible in Figures [5]^a) and (c), in that the scaled width v j Av of the 
t = ridge at low Av is about the same of that at high Av. In Figure [6] we also notice a 
systematic difference between the ridge widths of the two clusters, a feature that is already 
visible from the collapsed echelle diagrams of Figures (5^b) and (d). 

To see whether the difference in ridge width between the two clusters and between stars 
with different Av arises from the difference in temperature of the stars as contemplated by 
Chaplin et al] ( |2009[ ); |B~audm et al.| ( |2011| ); |Appourchaux et aI] ( |2012| ); [Belkacem et al.] ( |2012[ ), 
we plot our measurements of FWHM as a function of T e g in a log-log scale in Figure [7} This 
shows indeed that all the cluster stars follow an almost common trend, which supports that 
the observed difference in ridge width largely follows the difference in temperature. We also 
show the linewidth measurements of a sample of main sequence (MS) and subgiant stars 



(blue diamonds) from Appourchaux et al. (2012), where we have taken temperatures from 



Casagrande et al. (2006, 2010). It is remarkable how well all the stars are aligned in Figure [7j 
Note that our measure of the ridge width only provides an upper limit to the 'true' mode 
linewidths because of the slight curvature of the ridges in the echelle diagram. The fit to the 
linewidths across all stars is represented by an exponential function 



/7k- 5777K\ 
6XP ( Tq ) ^ 



(13) 



where r = 1.39 ± 0.10/iHz and T = 601 ± 3K (dashed black line). A detailed study 
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using linewidths found by direct mode fitting (peak bagging) of MS and red giant stars (e.g. 



Chaplin et al. 2009 Baudin et al. 2011; Appourchaux et al. 2012) goes beyond the scope of 



this work. The power law fit with a background component proposed by Appourchaux et al. 



(2012) is here added for comparison in its range of validity (5300 K— 6800 K, dot-dashed blue 
line). However, we can conclude that our measurements, combined with Kepler results on 
MS and subgiant stars, follow a single exponential trend with temperature. 



5. Mixed modes 



Mixed modes have the great advantage of being sensitive to the core structure, while 
at the same time being observable at the surface. They were recently used as a way to 



successfully distinguish between RC and RGB stars (Bedding et al. 2011; Mosser et al. 



2011a). Although their amplitude is lower than of pure p modes (Mosser et al. 2012a), 



long datasets enable us to identify many of them due to their long lifetimes (Dupret et al. 



2009). Even in cluster red giants, which are generally fainter than the Kepler field stars, we 
can detect many mixed I = 1 modes in the best cases. The main features of mixed modes 



relevant for the analysis presented in this work are discussed in Section 5.1 while our results 



on their period spacings are described in Section 5.2 



5.1. Period spacings of mixed dipole modes 



While p modes are equally spaced in frequency, pure g modes are approximately equally 



spaced in period, following the asymptotic relation (Tassoul 1980; Christensen-Dalsgaard 

Alio 



2011) 



n 



nt 



In + a) 



(14) 



Here, n and £ are the radial order and the angular degree of the mode, a is a small constant 
and An is the period spacing, expressed as 

-l 



An n = 2n' 



N 



dr 



(15) 



The integral is taken over the cavity in which the g modes propagate, and N is the Brunt- 



Vaisala frequency. From Equation (14), the period spacing of dipole g modes is given by 



AP g = AU /V2, 



(16) 
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which appears the most interesting quantity to investigate because of the strong coupling 
between p and g modes for t = 1 (iDupret et all 120091 iMontalban et al.||2010| IStellol 12011 



Bedding 2011). 



However, in contrast to the large separation for p modes, the period spacing of pure 
g modes, AP g , cannot always be directly measured in red giants because all the non-radial 



modes are mixed in the red giant phase ( Christensen-Dalsgaard 2011). Nevertheless, from 



recent studies it seems to be possible to infer n in some cases (see Bedding et al. 2011 



Mosser et al. 2012b). Fortunately we can readily measure the period spacing of the mixed 
modes, AP Q b s , which can serve as a proxy for AP g . AP Q b s is lower than AP g by about a 



factor of 0.6-0.8 (e.g. see Bedding et al. 2011 Mosser et al. 2012b) 



5.2. AP bs-Az/ diagram 



As mentioned in Sections [T] and [3j AARG measures period spacings using the approach 



used by Bedding et al. (2011). As a first step, it modifies the power spectrum for each star 
by erasing the regions containing all the £ = 0,2 modes, whose positions come directly from 
our analysis of p modes (see Section [3]). This new power spectrum shows only i = 1 mixed 
modes, and possibly some low amplitude £ = 3 modes, and is then expressed in period 
rather than frequency. The power spectrum of this power spectrum is then calculated, 
which is converted back into period. To obtain a first guess for the observed period spacing, 
a manual marking of the position of the excess of power is required. Finally, a Gaussian fit 
to the selected hump of power provides our measure of AP obs . 

We were able to measure AP Q b s for a total of 53 stars in our sample: 27 from NGC 6791, 
4 from NGC 6811, and 22 from NGC 6819. The results are shown in Figure [8j where RC 
and RGB stars are marked with open and filled symbols, respectively. The dashed black line 
is the limit set by the frequency resolution. Overlaid are theoretical tracks for 1.2 M e and 
1.6 Mq at near-solar metallicity (Z = 0.017), as calculated by White et al. (2011a), which 
are representative of the RGB stars of NGC 6791 and NGC 6819, respectively. Using stellar 
models, we verified that changing the metallicity over the range spanned by the two clusters 
has no significant effect on AP g for RGB stars. 

The tick marks at the top (RC) and bottom (RGB) represent stars for which the period 



spacing could not be clearly measured by our analysis, classified by Stello et al. (2011b) using 



the CMD. We see that the fraction of stars with measured period spacings is much higher 
for RC stars than for RGB stars even after taking into account the limit set by the frequency 
resolution. In particular, for NGC 6791 these fractions are ~36% (RGB) and ~73% (RC), 
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Fig. 8. — Period spacings of the three clusters NGC 6791 (red circles), NGC 6811 (blue 
diamonds), and NGC 6819 (green squares). Open symbols represent RC stars, while filled 
symbols are RGB stars. Tracks for 1.2 M Q (thick red line) and 1.6 M (thick green line), and 
Z = 0.017, are shown (White et al. 2011a[ ). Tick marks drawn at the top (RC) and bottom 
(RGB) refer to stars that could not be identified with our period spacing analysis, and are 
colored according to the notation adopted in the rest of the paper. The dashed black line 
represents the minimum period spacing one can measure with a 19 months-long time-series. 
Special cases discussed in Section 52 and listed in Table [T] are marked. 
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while for NGC 6819 they are -10% (RGB) and ~86% (RC). This could be explained by a 



weaker coupling between the p-mode and g-mode cavities for the RGB stars (Dupret et al. 



2009), which makes the resonances narrower in frequency, resulting in a smaller number of 



observable mixed modes. 



5.3. Discussion of special cases 



The stars labeled from A to D (KIC 5112361, KIC 4937770, KIC 5024414 and KIC 
5024476) are outliers in the AP ohs -Au diagram (Figure @), while KIC 2437103, KIC 2437589, 
and KIC 5024404 have period spacings that imply a different stage of evolution to the one 
based on the CMDs flStello et al.|2011b| ). KIC 9716522 represents a star on its way up towards 



the AGB, as already noted by Stello et al. (2011b) and now supported by our measurement 



of its high period spacing in agreement with that of other He-burning stars (Figure [8]). 

All these highlighted stars are also marked in Figures [3] and |4j and in the CMDs of 
Figure [9] (as derived by Stello et aL]|2011b" ). We also list these stars in Table [TJ together with 
all their asteroseismic parameters derived in this work. To further support the discussion 
presented below, we derived the masses of all stars near the RC in the CMD, including the 
outliers A— D and the stars that we have marked as 'likely evolved RC (Figure [8]), which 
are also shown in Figure [9] and listed in Table [TJ To estimate the masses we use the scaling 
relation 

M 1 \ [ Av \ f T eS 



where we adopted v. 



max,© 



Mq \ Z/ max ,0 

= 3100 fjHz, Av Q 



3/2 



AU Q J V T eff,0 

135 /iHz and T e ff i0 



(17) 



5777 K (e.g. Miglio et al. 



2012). The result is shown in Figure 10, with masses plotted against V magnitude and 1-er 



error bars overlaid. Blue lines represent the mean masses of RC stars (solid) and their 1-er 
uncertainties (dashed), as derived by Miglio et al. (2012) by adopting Equation (17). To 



provide corrected estimates of mass for clump stars, the Av scaling relation was corrected 
by 2.7 % and 1.9 % for NGC 6791 and NGC 6819, respectively, according to the study by 



Miglio et al. (2012). 



We first discuss possible causes for the outliers, A— D. All four are potentially binary 
stars. Three of them (A, C, and D) are listed as binary stars in the radial velocity study 



by Hole et al. (2009), and the fourth star (B) shows a low oscillation amplitude, which 



could be indicative of a binary star, as argued by Stello et al. (2011a). All four stars also 



appear relatively blue in the CMD (Figure [9|. Stars A and B fall below the RC in the 
CMD and are in line with the rest of the RGB stars in Figures [3] and [4], suggesting that 
they are RGB stars with no clear sign of an abnormal mass (Figure |4]( a) and (b)). This is 
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confirmed by our estimate of their masses according to Equation (17) (Figure 10), whose 



values are similar to the average mass of the RGB stars of NGC 6819 found by Basu et al. 



(2011) (see Section 4.2). Binarity seems like the most plausible explanation for their B — V 
colors being lower than the other RGB stars. However, their power spectra do not show 
oscillations from two components, and their higher-than-expected AP obs is therefore difficult 
to explain. Perhaps it could come from a different core structure of these stars caused by 
binary interaction. We note that the stars do not seem to be the result of a merger event, 
given their apparently 'normal' masses. In conclusion, stars A and B are most likely both 
binaries, with one component on the RGB, whose seismic signal we detect, and a fainter 
less-evolved component. Stars C and D have luminosities typical to that of the RC. Our 
measurement of AP obs suggests that the stars indeed belong to the RC. The position of 
the stars in the sequence of He-core burning stars going from low mass (low Au) to high 
mass (high Az/) spanned by the three clusters indicates that stars C and D have higher 
masses than the other RC stars in NGC 6819 (Figure |8]). This is confirmed by our estimate 



of their masses (Figure 10). Our conclusion that they are high- mass RC stars is in good 



agreement with Rosvick & Vandenberg (1998), who mention these stars along with others 



with this position in the CMD to be potential descendants of blue stragglers, meaning that 
they experienced mass transfer and therefore have a component with a mass significantly 
higher than the cluster's turn-off mass. 

Concerning the next three stars, our period spacing analysis shows that KIC 2437103 
(AP obs = 306 s) is an RC star, and KIC 2437589 (AP obs = 39 s) is an RGB star, as argued 

is an RC star. It seems that 



by Miglio et al. (2012), and KIC 5024404 (AP ( 



obs 



182 s 



KIC 2437589 is an evolved blue straggler in the RGB phase, as suggested by Brogaard et al. 



(2012). This would explain its unusual position in the CMD (top panel of Figure^]), and is 



also supported by a mass of about 1.7 M , as derived from Equation (17), greater than the 
mass of the other RGB stars of the cluster. 

Lastly, six stars (two in NGC 6791 and four in NGC 6819), were found to be possible 
candidates for RC stars that are starting to evolve towards the AGB. We list them as 'likely 
evolved RC in Table [TJ Our suggestion arises from our measurement of their AP obs , which 
corresponds to that of RC stars, and from their Az/, which is lower than that of the other 



RC stars. Their masses (Figure 10) are similar to that of the average RC star which, in 
combination with their lower Az/, confirm that they have a radius significantly larger than 
the other RC stars. 
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Fig. 9.— CMDs of the clusters NGC 6791 (top panel), NGC 6811 (middle panel), and NGC 
6819 (bottom panel) as derived by Stello et al. (2011b). Both RC and RGB stars are shown, 
with open and filled symbols respectively, according to the classification obtained by the 
membership study of Stello et al. (2011b) and our analysis of period spacings. Stars marked 
with labels represent special stars discussed in Section and listed in Table [T] Isochrones 
are shown for all the clusters (solid lines, see Stello et al. 2011b, for details). 
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Fig. 10.— Mass of stars near the RC of NGC 6791 (red circles) and NGC 6819 (green 
squares) with applied correction in the Av scaling of 2.7 % and 1.9 % respectively (see 



Miglio et al. 2012). Error bars show 1 — o uncertainties derived according to Equation rtl7 



Outlier stars A— D and 'likely evolved RC stars discussed in Section and listed in Table [T] 
are marked. Blue lines represent the corrected mean masses of RC stars (solid) and their 



1-cr uncertainties (dashed) derived by Miglio et al. (2012). 



6. Summary and Conclusion 



To summarize and conclude on the main results of our analysis: 



1. The fit of the e-Av relation to the RGB stars of our sample, computed using Equa- 
tion (|6]), is compatible with the result derived by Mosser et al. (2011b), although it 
deviates slightly towards higher values of Au, where our sample has more stars and 
benefits from longer observations. Our fit is almost indistinguishable from that by 



Kallinger et al. (2012), which was based on more than 900 field red giants observed for 



a similar length of time. Moreover, we tested a power-law form of the e-Au relation 
and found it to provide a very similar fit to that derived from Equation (|6]). Lastly, the 
average e of clump stars appears to be significantly different from that of their RGB 
counterparts for both NGC 6791 and NGC 6819, a result in agreement with previous 
findings on field RGs. 

2. The linear fits to the <5z/ 2-Az/ relation for the RGB stars of our sample, given by 
Equations and (J8j, appear to be compatible within a few percent with the results 
by Huber et al.| ( |2010 ) and Kallinger et al. (2012) on field red giants. A direct measure 
of the mass-dependence for the small spacings 5z/ 2 and dfoi is derived for the first time 
for cluster stars. The result indicates that Suqi is about three times more sensitive to a 
mass difference than 5vq2- The mass-dependence for 8vq 2 is compatible with the results 



by Kallinger et al. (2012) on field red giants. Furthermore, both Suq2 and Svqi show 



dependence on mass that is qualitatively in agreement with theoretical studies of red 



giant stars by Montalban et al. (2010). 



3. It is notable that the RC stars of NGC 6791 behave differently from those of NGC 6819 
for both <5z/ 02 and Su 01 , as visible in Figures [4] and |5j We quantified this unexpected fea- 
ture through the difference in {5v i) between RC and RGB stars, which is significantly 
different from one cluster to the other for both the small spacings. As discussed in Sec- 



tion |4.3[ further theoretical investigations concerning differences on mass, metallicity 
and rotation between the two cluster populations, are required for the full interpreta- 
tion of our results. 

4. The positions of the ridges in the collapsed ensemble echelle diagrams (Figure [5]) con- 
firm the results from Huber et al. (2010), with the £ = 1 and 2 ridges moving away from 
the £ = ridge as the stars evolve from the H-shell to the He-core burning phase. The 
position of the £ = 3 ridges, hence of the average small spacings 5u 03 (Figures |5](b) and 



(d)), is also in agreement with results by Bedding et al. (2010); Huber et al. (2010); 
Mosser et all pOlTbl; iKallinger et all d25l2l) on field red giants. The FWHM of £ = 
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ridge, which represents an upper limit of the mode linewidths, increases with Au for 



both NGC 6791 and NGC 6819, a result that agrees with the studies by Huber et al. 



(2010) and by Kallinger et al. (2012) on field red giants. A systematic difference of 



the FWHM between the two clusters is shown, which is largely explained by the tem- 
perature dependence of mode linewidths (Figure [7]), as discussed in Section 4.5 This 
result shows the first evidence for an exponential correlation between mode linewidth 
and temperature in red giants, which is consistent with extrapolating Kepler results 



for main-sequence and subgiant stars derived by Appourchaux et al. (2012) 



5. The analysis of period spacings using the method described by Bedding et al. (2011) 



allowed for the successful identification of almost half of the stars in our sample as 
either H-shell or He-core burning red giants. The fraction of RGB stars with clearly 
detectable period spacings is much lower than for RC stars, as seen in Figure |8j and 
could be due to a weaker coupling of the p- and g-mode cavities in RGB stars. We 
see a number of outliers in Figure [8] which require further investigations. It is likely 
that all of them are binaries and two of them appear to be evolved blue stragglers, 



as suggested by Rosvick & Vandenberg (1998) and supported by their higher masses 



(Figure 10). In addition, our analysis confirms the suggestion by Miglio et al. (2012) 
that stars KIC 2437103 (AP obs = 306 s) and KIC 2437589 (AP obs = 39 s), are an RC 
and an RGB star, respectively, and that KIC 5024404 (AP Q b s = 182 s) is an RC star. 
Lastly, we find a number of possible candidates for evolved RC stars in both NGC 6791 



and NGC 6819, as suggested by our measurement of their masses (Figure 10), which 



indicate they have a radius larger than the other RC stars. The special cases discussed 



in Section 5.3 represent potentially interesting targets for detailed theoretical modeling. 
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